Nonlinear dynamics of electromagnetic pulses in cold relativistic plasmas 
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In the present analysis we study the self consistent propagation of nonlinear electromagnetic pulses 
in a one dimensional relativistic electron-ion plasma, from the perspective of nonlinear dynamics. 
We show how a series of Hamiltonian bifurcations give rise to the electric fields which are of relevance 
in the subject of particle acceleration. Connections between these bifurcated solutions and results 
of earlier analysis are made. 



I. INTRODUCTION 

Propagation of intense electromagnetic pulses in plas- 
mas is a subject of current interest in a variety of ar- 
eas that make use of the available modern laser tech- 
nologies, among which we include particle and photon 
acceleration, nonlinear optics, laser fusion, and others 
IB IE [J @- Intense electromagnetic pulses dis- 
place plasma electrons and creates a resulting ambipo- 
lar electric field with the associated density fields. Un- 
der appropriate conditions all fields act coherently and 
the pulse keeps it shape. Studies on pulse localization 
have been performed in a variety of forms to unravel 
the corresponding numerical and analytical properties. 
Kozlov et al. Q investigate numerically propagation of 
coupled electromagnetic and electrostatic modes in cold 
relativistic electron-ion plasmas to conclude that small 
and large amplitude localized solutions can be present. 
Mofiz & de Angelis [8| apply analytical approximations 
to the same model and suggest where and how localized 
solutions can be obtained. Ensuing, more recent works 
provide even deeper understanding as various features 
are investigated, like influence of ion motion in slow, ion 
accelerating solitons [jj , existence of mo ving so litons [U , 
existence of trails lagging isolated pulses [lfllU and oth- 
ers. Some key points however remain not quite under- 
stood, like the way small amplitude localized solutions 
are destroyed; when isolated pulses are actually free of 
smaller amplitude trails; and more specific properties of 
the spectrum of stronger amplitude solutions, to mention 
some. Those are issues of relevance if one wishes to es- 
tablish the existence range and stability properties of the 
localized modes. 

In the present paper we shall turn our attention to 
small amplitude solitons propagating in undcrdense rar- 
ified plasmas, since these kind of solitons may be of rele- 
vance for wakefield schemes. In doing so we shall follow 
an alternative strategy, other than the direct integration 
of the governing equations which has been the standard 
approach so far. We intend to examine the problem with 
techniques of nonlinear dynamics [T^ | . A canonical repre- 
sentation shall be constructed in association with several 
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tools of nonlinear dynamics like Poincare maps and sta- 
bility matrices. This strategy naturally provides a clear 
way to investigate the system since we intend to estab- 
lish connection between the pulses of radiation and fixed 
points of the corresponding nonlinear dynamical system 
|l3j | . Several facts are known and we state some which are 
of direct relevance for our analysis: small amplitude soli- 
tons are created as the wave system becomes modulation- 
ally unstable at an upper limit of the carrier frequency 
and cease to exist beyond a lower limit of this carrier 
frequency. Not much is known on how solitons are de- 
stroyed at the lower boundary and we examine this point 
to show that a series of nonlinear resonances and bifurca- 
tions are responsible for process. A related relevant prob- 
lem is when isolated pulses are actually free of smaller 
amplitude trails and this has to do with the existence of 
wakefields following the leading wave front which is of 
relevance for particle acceleration, for instance. Those 
are basic issues if one wishes to operate the wave system 
under conditions suited for particle acceleration, and our 
purpose with the present paper is to contribute towards 
the analysis of these aspects. 



II. THE MODEL 

We follow previous works and model our system as 
consisting of two cold relativistic fluids: one electronic, 
the other ionic. Electromagnetic radiation propagates 
along the z axis of our coordinate system and we rep- 
resent the relevant fields in the dimensionless forms 
eA(z,t)/m e c 2 — > A(z,i) for the laser vector potential, 
and e4>{z,t)/m e c 2 — > <f){z,t) for the electric potential. 
— e is the electron charge, m e its mass, and c is the speed 
of light; mi will denote ionic mass when appropriate. In 
addition, we suppose stationary modulations of a circu- 
larly polarized carrier wave for the vector potential in the 
form A(z,t) — il>(0 [ x s'm(kz — cot) + y cos(fc;z — cut)] with 
£ = z — Vt, whereupon introducing the expression for 
the vector potential into the governing Maxwell's equa- 
tion one readily obtains V = c 2 k/ui. V could be thus read 
as a nonlinear group velocity since we shall be working in 
regimes where oj and k are related by a nonlinear disper- 
sion relation. Manipulation of the governing equations 
finally takes us to the point where two coupled equations 
must be integrated - one controlling the vector potential, 
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and the other the electric potential @, : 
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where the primes denote derivatives with respect to £ = 

( tte/c)|, r e (0,V>) EE A /(l + ^)2-p(l + ^2) ) r .(^) = 

•/(r^T^p - T^V 5 ), ?7 = u 2 e /uj 2 , A* = me/mi, 
Vo ee V/c, and p e 1 - Vq 2 , with cj 2 ee 47m e e 2 /m e as 
the plasma frequency, and Tig — ^"i 

the equilibrium 

densities. We further rescale w/ck — > w and u e /ck — > oj c 
in Vo, ?7 and p, which helps to simplify the coming investi- 
gation: ?7 preserves its form, Vq — > 1/w, andp — > 1 — 1/lu 2 . 
A noticeable feature of the system Q - (J2J is that it can 
be written as a Hamiltonian system of a quasi-particle 
with two-degrees-of-frecdom. Indeed, if one introduces 
the momenta P^ = ip' and ee —<fi'/p, the equations 
for ip and takes the form 



ip' = dH/dP^, = -dH/dijj, 
0' = dH/dPj,, P; = -dH/dh 



where the Hamiltonian H reads 
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H is constant since it does not depend on the "time" 
variable £. Its constant value, let us call it E, can 
be calculated as soon as the appropriate initial condi- 
tions are specified. In our case we shall be interested 
in the propagation of pulses vanishing for |£| — > oo, so 
we know that conditions P^ = Pq = = = must 
pertain to the relevant dynamics, from which one con- 
cludes that E — (Vo/p) 2 + Additional conditions 
arise from the presence of square roots in the Hamilto- 
nian; the dynamics lies within regions where simultane- 
ously r 2 , r 2 > 0. Combining these inequalities with the 
boundary conditions one is led to conclude that the entire 
dynamics must evolve within the physical region 
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if p > 0. If p < there is no restriction, but we shall see 
that only positive values of p are of interest here. We can 
also evaluate the linear frequencies of laser and wakefield 
small fluctuations in the form 
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where 



-l/r) + l/p(l + (x), and ft 2 ee (1 + ^/V 2 . (8) 



The potential oscillates with a real frequency f2^ which 
can be shown to convert into u) e (l + A*) 1 / 2 if dimensional 



variables are used for space and time. As for the vector 
potential, to reach high-intensity fields from noise level 
radiation, instability must be present, which demands 
fi 2 > and, consequently from relation JSJ, 



1< lo 2 < 1 +uj 2 e (1 + //), 



(9) 



so p > 0. 



The threshold = can be rewritten in the form lu = 



ee y/l + uj 2 (1 + /i), where is the linear dispersion 
relation for electromagnetic waves. What we expect to 
see are small amplitude waves when u> is slightly smaller 
than w*, with amplitudes increasing as we move farther 
from the threshold. In addition to that, another feature 
worth of notice must be commented. If one sits very close 
to the threshold, amplitude modulations of the laser field 
are tremendously slow, while the oscillatory frequency 
of the electric potential remains relatively high. The 
resulting frequency disparity provides the conditions for 
a slow adiabatic dynamics where, given a slowly varying 
-0, always accommodates itself close to the minimum 
of 

1/(0, 0) = -v /p 2 [r e (0, 0) + »- x n{4>, VO] , (10) 

the "minus" sign on the right hand side accounting for 
the negative effective mass of as seen in Eq. (JSJ; note 
that m m of Eq. © refers to the smallest available and 
not to the minimum of U . When ip = 0, a condition to 
be used shortly in our Poincare plots, U has a minimum 
at = which is thus a stable point in the adiabatic 
regime. As one moves away from the threshold, faster 
modulations and higher amplitudes may be expected to 
introduce considerable amounts of nonintegrable behav- 
ior and chaos into the system. This kind of perspec- 
tive agrees well with the result of previous works where 
adiabatic regions have been interpreted to be essentially 
associated with small amplitude quasineutral dynamics 
pj . One of our interests here is to precisely see how the 
adiabatic dynamics is broken as one moves deeper into 
nonintegrable regimes. An additional fact must be ob- 
served as one searches for adiabatic solutions and this 
has to do with how close to the minimum of U on must 
sit so as to find these adiabatic solutions. The corre- 
sponding discussion parallels that on wave breaking of 
relativistic eletrostatic waves. First of all note that if 
we do not set right at the respective minimum of U, 
the electric potential will oscillate around the minimum 
which will be itself displaced due to the action of the 
slowly varying ip. Again when ip = 0, inequality (JBJ re- 
veals that must lie in the range mm = ^fp — 1 < to 
4>max — 1/^(1 — \/p) > 0- Not all these values are how- 
ever actually allowed in adiabatic dynamics. Oscillations 
will occur consistently only if the orbit is free to wander 
to the right and left hand sides of the minimum = 
and this can only happen when the oscillating orbit is 
entirely trapped within the attracting well of U. U < 0, 
and a quick calculation shows that 

U(4> min Y - U((f> max ) 2 = 2VP(1 - A< 2 )x 
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FIG. 1: Oscillating (7) and wave breaking (77) regions for the 
electric potential at ip — 0. AU is denned in the text. 
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so U(4> rnax ) > U((f> m i n ), which sets a limit to cyclic or- 
bits: <p must be such that the corresponding potential 
will never be above the level U(4> m in)- To illustrate all 
these comments, the reader is referred to Fig. ^ where 
the potential AU = U((p,tp = 0) - U(<j> = 0, ip = 0) is 
represented: orbits of region I, <p m i n < <j> < 4>, will os- 
cillate back and forth, but orbits in region II eventually 
reach 4>min where r e — > 0. 

Since it can be shown that the electronic density de- 
pends on r e in the form n e ~ r~ x 0, @|j break down 
of the theory indicates wave breaking on electrons. Also 
shown in the figure is the wave breaking energy 
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separating regions I and 77. Our conclusion is that even 
with extremely slow modulations, oscillations of (p must 
be limited so as to satisfy the conditions discussed above. 
Not only that, but the very same figure suggests how non- 
integrability affects localization of our solutions: as one 
moves away from adiabaticity and into chaotic regimes, 
trajectories initially trapped by U may be expected to 
chaotically diffuse towards upper levels of this effective 
potential, escaping from the trapping region, approach- 
ing E w b r and eventually hitting the boundary at (p m in 
or, in general, attaining r e — for ip ^ 0. If this is so, 
we have an explanation on how small amplitude solitons 
are destroyed, one of the issues of interest in the subject 
[B| . We now look at the problem with help of methods of 
nonlinear dynamics. 



III. ANALYSIS WITH NONLINEAR 
DYNAMICS 

We introduce our Hamiltonian phase space in the 
form of a Ponicare surface of section mapping where the 
pair of variables (<j>, P^) is recorded each time the plane 
ip = is punctured with P^ < 0. Once we have de- 
fined the map this way, we can also investigate the ex- 
istence and stability of periodic solutions of our coupled 
set of equations with the aid of a Newton-Raphson al- 
gorithm. The Newton-Raphson method locates periodic 
orbits and evaluates the corresponding stability index a 
which satisfies \a\ < (>)1 for stable (unstable) trajecto- 
ries |l4j . Parameters are represented in a form already 
used in earlier investigations on the subject: we first 
set a numerical value for Vq and then obtain to = 1/Vo 



which must be larger than the unity as demanded by 
condition 0. However, we shall keep V close to the 
unit, and thus lo slightly larger than one, so as to rep- 
resent wave modes propagating nearly at the speed of 
light. This is the convenient setting if one is interested 
in fast electron acceleration by wakefields. After V is 
established, the electron plasma frequency is calculated 
as = rjuj 2 , r\ satisfying condition (J5J again. We note 
that 77 = uj 2 Jlj 2 = V 2 u 2 e = V 2 ul nonscaled /c 2 k 2 , so hold- 
ing Vo constant while increasing 77, is entirely equivalent 
to the more usual practice of holding Vo and the orig- 
inal U! ei while decreasing k and the original to. In all 
cases analyzed here we take fi = 0.0005 as in Kozlov et 
al 0- In addition to that, we shall take V Q = 0.99 to rep- 
resent the high speed conditions of wakefield schemes. 
Now a crucial step is this: since isolated pulses cannot 
be seen in periodic plots we alter slightly the energy E 
to E = V /p 2 (l + 1/ fi) (1 + e), e < 1 so the vanishing 
tail P,/, = Pfj, = ip = <f> = is avoided. With this maneu- 
ver we convert isolated pulses into trains of guasi-isolated 
pulses, a situation amenable to the use of nonlinear dy- 
namics and the associated periodic plots; periodicity is 
in fact physically meaningful if pulses result from peri- 
odic self-modulations of initially uniform modes |l5|. In 
all cases we make sure that as e — * the trains go into 
individual packets and convergence is attained. The in- 
stability threshold for the vector potential is obtained in 
the form 77* = p/(l +/i) = 0.0198 so u> v <C oj as it must be 
in the underdense plasmas. To investigate the adiabatic 
regime of the relevant nonlinear dynamics we examine 
phase portraits for 77 slightly larger than 77* . In panel (a) 
of Fig. Ewe set 77 = 1.00001 77*. 

such a relatively small departure from marginal sta- 
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FIG. 2: (a) Phase plot near the modulational instability 
threshold, with r\ = l.OOOOlr;*; (b) stability index versus r/; 
(c) phase plot after the inverse tangency seen in panel (b), 
with T) = 1.0001??,. e = 10" 11 . 



bility, modulations are slow with 3> adiabatic 
approximations are thus fully operative and what we see 
in phase space is just a set of concentric KAM surfaces 
rendering the system nearly integrable. The central fixed 
point corresponds to an isolated periodic orbit since it 
represent a phase locked solution that return periodi- 
cally to if) — <fi — > 0, and the surrounding curves de- 
pict regimes of quasiperiodic, nonvanishing fluctuations 
of 4>. Resonant islands are already present but still do 
not affect the central region of the phase plot where the 
solitary solution resides. In general we have observed 
that increasingly large resonance islands are present away 
from the central region. When 77 grows the behavior of 
the central fixed point can be observed in terms of its 
stability index represented in Fig. E^b). The index os- 
cillates within the stable range initially, which marks the 
existence of a central elliptic point near the origin. The 
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FIG. 3: Dynamics as represented in the versus $ space: 
(a) r\ = I.OOOOIt?,; (b) 77 = 1.000217?*. e wbr = X eE wbr /( X e + 

^wbr ) - 



stability index however finally reaches a = +1 as indi- 
cated in the figure, beyond which point no orbit is found. 
This indicates a tangent bifurcation with a neighbouring 
orbit which terminates the existence of the central point 
[lfil | . Immediately after tangency, the phase plot at ip = 
is still constricted to small values of </> as seen in Fig. ^c) 
where rj = 1.000177*. Larger values of rj cause diffusion 
towards upper levels of U (0) and we can see that in Fig. 
(3), where we investigate the behaviour of the energy 



(13) 



corresponding to the electrostatic field (p. Instead of 
working directly with the form (|13fl we represent diffusion 
in terms of compact variables 



: X4- + W 



(14) 
(15) 



where Xe and \4> represent the scale above which the 
corresponding variables are compactified. 

This kind of choice allows us to represent in the same 
plot the very extensive variations of energy and electric 
potential, without deforming these quantities when they 
are small, near their initial conditions. We found it con- 
venient to use Xe — X<t> — 0.0001 to discuss diffusion. 
In Fig. |3{a) we take rj — I.OOOOI77* so we are in the 
regular regime; as expected, no diffusion is observed and 
the quasi-particle stays near its initial condition = 0, 
cf> = 10 -8 . For 77 = 1.0002177* as in panel (b), the central 
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FIG. 4: "Time" series for the electric field \dcf>/d£\ for r\ — 
1.0004t?«. 



As rj increases past the inverse tangent bifurcation but 
prior to full destruction of isolating KAM surfaces, one 
reaches a regime of periodical returns to ip = 0, although 
in the presence of a slightly chaotic 4> motion. Those 
cases where ip = but 4> ^ 0, correspond to quasineutral 
ip pulses accompanied by trails of </> activity as described 
in Refs. 0] and 0] • We see that trails can be regular or 
chaotic. Finally, for large enough rj's, KAM surfaces no 
longer arrest diffusion and wave breaking does occur as 
r e — > 0, as we have checked. At this point adiabatic mo- 
tion is lost and this is likely to correspond to that point 
where small amplitude solitary solutions are entirely de- 
stroyed as commented in Refs. [j| and 0. 



IV. FINAL CONCLUSIONS 



fixed point no longer exist. In addition to that, KAM 
surfaces no longer isolate the central region of the phase 
plot and diffusion is observed. The quasi-particle moves 
toward E w \, r and eventually arrives at this critical energy 
producing wave breaking on electrons. At this point the 
simulation stops with the electron density diverging to 
infinity. Diffusion is initially slow and becomes faster as 
energy increases. One sees voids in the diffusion plots 
which correspond to resonant islands in the phase space, 
so as diffusion proceeds the quasi-particle escalates along 
the contours of the resonances that become progressively 
larger as already mentioned - this is why the process is 
initially slow, becoming faster in the final stages. For 
larger values of r\ no resonance is present and the quasi- 
particle moves quickly toward E w f, r . In case of panel (b) 
one can still see various pulses before wave breaking, but 
when rj is so large that resonances are no longer present, 
wave breaking can be instantaneous. We finally note the 
following relevant fact. For V — > 1, it is known that 
the amplitude of the electromagnetic pulses are small [f| . 
But as one goes beyond the adiabatic regime, our dis- 
cussion on diffusion allows to conclude that even small 
initial pulses eventually reach very high values, which 
provides the condition for formation of strong electric 
fields with the corresponding implications on particle ac- 
celeration. We illustrate the feature with a final figure, 
Fig. 01 where, in a diffusive regime with ?y = 1.0004^„, 
the electric field —<f>' = —pP<p is shown to evolve from 
small values near initial conditions to the limiting wave 
breakin g value which agrees with the calculated value 

\<//\ ~ y/2UJU~ e ~ 3.5. 

We read all these features as it follows. For small 
enough 77's there are locked solutions representing iso- 
lated pulses coexisting with surrounding quasiperiodic 
solutions where <f> does not quite vanish when ip does. 



To summarize, we have used tools of nonlinear dy- 
namics to examine the problem of wave propagation in 
relativistic cold plasmas, discussing underdense regimes 
appropriate to wakefield schemes. Nonlinear dynamics 
provides a unified view on the problem, thus allowing to 
address simultaneously several relevant questions. In this 
paper we have kept our interest focused on weakly non- 
linear modes where a transition from adiabatic to nonin- 
tegrable dynamics was observed. Starting with very low 
amplitude regimes near the onset of modulational insta- 
bility, one has either isolated pulses or pulses coexisting 
with regular </> trails. As one increases rj, thus moving 
away from the onset, pulses with slightly larger amplitude 
exist but are never fully isolated since tangent bifurca- 
tions annihilate the central fixed point and create ubiq- 
uitous chaotic electrostatic trails. However, electrostatic 
activity is still surrounded by KAM surfaces and there- 
fore confined to small amplitudes. Now as one pushes 
amplitudes a little higher, isolating KAM surfaces are 
destroyed, pulses are no longer possible at all and wave 
breaking does occur. There are therefore three clearly 
identified regimes in the problem: (i) regular or adiabatic 
regimes where the dynamics is approximately integrable, 
(ii) a weakly chaotic regimes where chaos is present but 
chaotic diffusion is still absent due to the presence of lin- 
gering isolating KAM surfaces, and finally (iii) diffusive 
chaotic regimes where isolating KAM surfaces are absent. 
Thermal effects should be added All those issues are of 
significant importance for plasma accelerators and shall 
be developed in future publications. 
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